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Engineering semiconductors to enhance carrier multiplication (CM) could lead to increased pho- 
tovoltaic cell performance and a significant widening of the materials range suitable for future solar 
technologies. Semiconductor nanocrystals (NCs) have been proposed as a favourable structure for 
CM enhancement, and recent measurements by transient absorption have shown evidence for highly 
efficient CM in lead chalcogenide and CdSe NCs. We report here an assessment of CM yields in 
CdSe and CdTe NCs by a quantitative analysis of biexciton and exciton signatures in transient pho- 
toluminescence decays. Although the technique is particularly sensitive due to enhanced biexciton 
radiative rates relative to the exciton, k r gx > 2fc™ d , we find no evidence for CM in CdSe and CdTe 
NCs up to photon energies hu) > 3E g , well above previously reported relative energy thresholds. 



Carrier multiplication (CM) in the form of impact 
ionization is a w ell-understood phenomenon in bulk 
semiconductors^!. The process, consisting of inelastic 
scattering of energetic charge carriers and valence elec- 
trons to create additional e-h pairs, normally has high 
energy thresholds and low efficiency due to momentum 
conservation requirements and competition from ultra- 
fast intraband relaxation. While this conventional bulk 
CM has found a particular application in avalanche pho- 
todiodes for single photon detection, efficient CM fol- 
lowing optical excitation could have a significantly wider 
impact in the area of solar energy conversion. 

In a typical photovoltaic cell with a single active layer, 
photon energy in excess of the bandgap is lost by rapid 
thermalization. The CM process, if efficient, could har- 
vest this excess energy into additional e-h pairs, boost- 
ing the maximum theoretical power conversion efficiency 
from 32% to >40%P^, and, more importantly, widen- 
ing the range of candidate materials for new solar tech- 
nologies to include previously ignored narrow-gap semi- 
conductors. Strongly confined semiconductor nanocrys- 
tals (NCs) have been proposed as candidate structures 
for efficient CM^ because of an anticipated relaxation 
of momentum conservation constraints^ and potential 
slowing of competing phonon-mediated intraband cool- 
ing due to the discrete electronic structure ("phonon 
bottleneck" )0. 

Recently, transient absorption (TA) measurements 
have shown evidence of efficient ultrafast CM in IR- 
emitting PbSe, PbS, and PbTe NCs above a photon 
energy threshold hu = 3-Ej 8 * 9 * 10 ! In extensions of the 
work, subsequent TA measurements have indicated that 
a single-photon could generate up to 7 e-h pairs in PbSe^ 
and that CM is similarly efficient in visible-emitting CdSe 
NCs above a Huj = 2.5E g thresholcP 1 . 

The conclusions of the TA measurements suggest new 
and unique underlying physics as well as some interesting 
questions. First-principles theories explaining the bal- 
ance of Coulomb coupling and phonon relaxation rates 
implied by the experiments have yet to emerge. At the 
same time, studies on intraband relaxation in CdSe and 



PbSe NCs at room temperature have found fast cooling 
dynamics that do not appear consistent with a phonon 
bottleneckPSEIl. r n addition, some aspects of the ex- 
perimental data are intriguing, such as similar CM ef- 
fects seen in PbSe and CdSe despite the very different 
state structures at threshold, and the observed linear de- 
pendence of CM yields on excess energy^. 

We present here an assessment of CM yields in CdSe 
and CdTe nanocrystals by transient photoluminescence 
(tPL). While complementary to TA in some ways, tPL 
is a background-free measurement better suited to the 
low excitation fluences necessary in CM studies. It is 
also more selective since it relates to the number of e- 
h pairs instead of single-particle state filling. Both TA 
and tPL capitalize o n the un ique, fast dynamics of mul- 
tiexciton (MX) state d 16 | 17 | ls l to isolate and quantify MX 
populations, and though more complicated to interpret, 
tPL becomes a useful and particularly sensitive technique 
when carefully calibrated. A recent study has analysed 
CM yields in CdSe NCs by tPL and concludes that there 
is agreement with previous TA determinations^. Their 
analysis and experiment differs from ours in several sig- 
nificant ways. We find instead that CM efficiency in CdSe 
and CdTe NCs is close to zero even for photon energies 
up to 3.lE g . Implications are discussed. 

Tunable UV excitation pulses were generated by non- 
linear mixing of the visible output and 3.1 eV remnant 
of an optical parametric amplifier (Coherent OPA 9400) 
pumped by a 250 kHz amplified Ti:sapph laser (Coher- 
ent RegA 9000). In a seperate experiment, the Ti:sapph 
was tuned to the red and 5.9 eV pulses were obtained by 
doubling its second harmonic. Room temperature hexane 
dispersions of NCs in lmm-path length cuvettes were ex- 
cited at 45° incidence. Emission was collected front-face 
and spectrally dispersed onto a streak camera (Hama- 
matsu C5680). We studied both organic ligand-capped 
CdSe or CdTe (core) and CdSe/ZnS or CdSe/ZnCdS 
(core/shell) overcoated particles^. 

Fig[l] shows PL decays of a representative sample of 
E g =2.0 eV CdSe NCs under weak and strong excitation 
at 3.1 eV and 5.6 eV (E g determined from the lowest 
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FIG. 1: (color online), (a) Absorption spectrum of a typical 
CdSe NC core sample used. Excitation energies employed in 
tPL are indicated by the arrows (inset) Detail of band edge 
absorption structure and emission spectrum, (b) Band edge 
luminescence dynamics of the sample at the two indicated 
excitation energies for peak pulse fiuence as noted. 



absorption feature). At high fiuence both decays show 
an additional fast component consistent with biexciton 
(BX) emission. Remarkably, and unlike data from RefHUI 
which shows a fast component under UV excitation, we 
find that the two low-fluence decays follow each other 
closely, suggesting that CM is less efficient than previ- 
ously reported. A quantitative determination of the CM 
yield requires first a careful characterization of the BX 
tPL signature, which we describe in detail. 

Results of an excitation- fiuence series at 3.1 eV are 
shown in Fig [2j Decays at low fiuence are dominated 
by exciton (X) emission. In most cases, X dynamics are 
multiexponential with contributions from NCs wit h dif- 
ferent non-radiative relaxation and trapping rate d 20 * 2 ^. 
The decays are adequately described by a biexponential, 
f(t) = cxp(-t/T Xs iow) + Cf ast exp(-t/Txfast), chosen 
for simplicity, where Txfast ~ 100-300 ps and Txsiow ~ 
1-10 ns. With increasing fiuence, early-time spectra show 
biexciton (BX) emission at the band edge and a fur- 
ther, blue-shifted feature corresponding to I P-IP emis- 
sion from higher multiexcitons (Figj2^i &[2]d) 1718 . The 
BX state then decays quickly with a size-dependant life- 
time tbx ~ 0.1-1 ns due to a fast non-radiative "Auger" - 
like coulomb process^. As shown in Fig|2]:, the mea- 
sured tPL decays are well described by a superposition 
of X dynamics and an additional single exponential BX 
component, a B x exp (—t/rBx) + axf(£p®. 

The relationship between the observed tPL decay am- 
plitudes, ax and clbXi and the underlying MX popu- 
lations was studied using a first-order kinetic model of 
MX relaxation, p' m (t) = -k m p m + k m+1 p m+1 , giving 
a x oc pZ l>0 k x d and 
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where p m is the relative population of NCs with m elec- 
trons and holes, p° n are initial values at t=0, k m are 
the MX decay rates, and and k B "x are the X and 
BX radiative rates. The asx term includes BX popu- 
lations formed by cascaded decay of higher MX states 




energy (eV) 



time (ns) 











740 uJ/cm 2 






250 fiJ/cm 2 






1 05 jiJ/cm 2 






21 uJ/cm 2 




time (ns) 

FIG. 2: (color online), (a) Transient PL spectra of the NCs 
in Fig.l collected from t — — 10 ps to t = 10 ps and (b) 
decays integrated from 1.89 to 2.02 eV under 3.1eV excita- 
tion with a peak fluences as noted, (c) A separate subset 
of our data, showing PL decays normalized at long times. 
The black lines are fits to the form asxe.~ t ^ TBX + axf(t) 
where tbx = 185 ps and f(t), the single X dynamics, are 
kept constanlpM] (d) (lines) Predicted initial NC popula- 
tions in an X or higher state (p m >o) or in a BX or higher 
state (p m >i) plotted against the peak average number of e-h 
pairs created, no = max{n(r)}. ( x >°) Scaled and fit ax and 
a,Bx PL components extracted from fits to measured decays. 



and is proportional to k B a ^ — fc™ d since BX emission is 
partially offset by correspondingly reduced X emission at 
early times^. The p° m are related to the incident laser 
power assuming poissonian photon absorption statistics 
and explicitly accounting for the inhomogcncous excita- 
tion profile of the beam: 



n(fy 



n(r) = j P (f)a 



where j p (r) is the measured photon flux at r, n(r) is the 
average number of absorbed photons per NC, and a is the 
absorption cross-section, treated here as an adjustable 
parameter. 

Band edge PL decays show growth and slow satura- 
tion of ax and asx that fit reasonably well to the ex- 
pected curves (see Fig[2]i)P^I From the fits we extract the 
IbxIo'X ratio expected at BX saturation (p° l=2 = p„ >0 ) 
and find a sample-dependant (asx I 'a>x)sat value in the 
range 3-6. This implies a substantially faster radiative 
rate of the BX relative to the X and leads to enhanced 
sensitivity of tPL for detection of small multiexciton pop- 
ulations, as is illustrated by the prominence of BX fea- 
tures in Figj2|;. We estimate k B a ^-/k r x d > 3-5, consistent 
with previous measurements on high quality NCs show- 
ing k r ££ /k r ^ d « 3^21. 

The observed faster radiative lifetime of the BX rela- 
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FIG. 3: (color online). Band edge PL decays of E g =2.07 eV 
(top) and E g =1.8Q eV (bottom) core/shell CdSe NCs under 
weak (no < 0.01) excitation at 3.1eV (black) and 5.6 eV (red). 
The 5.6 eV excitation corresponds to hu>=2.70E g and 2.95E g 
respectively, for which TA measurements at 6.2 eV predict 
CM yields of 22% and 50%^ 



tive to the X is explained by the electronic fine structure 
of these states. Because emission from the X ground 
state is spin-forbidden, X luminescence is relatively slow 
(£;™ d « 0.05 ns -1 ) and consists mostly of emissi on from 
thermally populated spin-allowed bright states^- 24 * 25 !. 
However, transitions from the BX ground state to some 
states in the X fine structure are predicted to be opti- 
cally allowed 25 , indicating that k r BX can be significantly 
larger than the 2k x ld value one would predict from a sim- 
ple carrier counting argument. 

The sensitivity of the tPL experiment was thus ex- 
ploited to estimate carrier multiplication yields in NC 
samples under UV excitation^. We find that, while sig- 
natures of multiexciton emission appear at high fluence, 
decays under weak 5.6 or 5.9eV excitation are close to 
indistinguishable from decays under weak 3.1eV (Fig. [3| 
even for large NCs where the excitation energy fno ex- 
ceeds 3E g . This contrasts with the measurements of Ref. 
1191 which show an additional fast component under UV 
excitation. Quantitative values shown in Fig|4]of the CM 
yield, the fraction of photo-excited NCs initially found in 
the BX state, y cm = p f* p J ^ , were obtained as the ratio 

y cm = {asx/ax) — fr° m fitting the UV-excited decays to 
the form a BX exp (-t/r BX ) + a x f(if^l. 

As seen in Fig [4] our results do not match the Hui/E g 
dependence of y cm found in TA measurements on CdSe 
NCs using 6.2 eV excitatioip] (FigQ. This could be 
because y cm is not only a function of the ratio tun/Eg, 
as suggested by TA on PbSe NOP, but also depends 
explicitly on huj. However, our interpretation of re- 
cently published tPL data at 6.2 eV^ using our esti- 



FIG. 4: (color online). Carrier multiplication yields un- 
der 5.6 and 5.9 eV excitation extracted from tPL decays, 
y cm = (agx X /ix) t ' Conservative (small) values {a B x / a x ) sat 
were used for each sample. 2er-wide error bars from repeated 
measurements are shown if larger than the symbols. The 
dashed line is the y cm vs. ftw/E g dependence found in the 
TA studies of RefQj] using hu> = 6.2eV. 



mated k r BX > ik r x d suggests a y cm of at most « 35% 
even at fko/E g w 3.2, instead of the 70% CM yield re- 
ported, which assumed fcg^ = 2k r ^^. tPL and TA 
assessments of CM yield in CdSe NCs therefore appear 
to disagree, and it is possible that CM in semiconductor 
NCs is generally less efficient and not as universal as has 
been thought. 

There are several reasons to expect a different CM en- 
hancement in II- VI semiconductors compared to the lead 
chalcogenides. The electronic state structure in wide gap 
NCs at energies > lE g in excess of the band edge is likely 
bulk-like with lcvcl-spacings that are small. On the other 
hand, the analogous electronic states in the lead chalco- 
genides might be more discrete in character because of 
the smaller E g and lighter effective masses. For these rea- 
sons, our experimental conclusions cannot be extended to 
the narrow band NCs, but the findings suggest the need 
to verify TA assessments of y cm in the lead chalcogenides 
with other techniques such as tPL. 

We note that the nature of a potential CM enhance- 
ment mechanism in NCs is itself not well understood. Al- 
though theoretical schemes have made progress towards 
explaining some of the reported phenomenology, these 
calculations have yet to quantitatively account for the 
balance of Coulomb and intraband relaxation processes 
leading to CM enhancement in a way that is clearly con- 
sistent with what is already known about carrier cooling 
and Auger rates in NCs. Calculations based on a tradi- 
tional impact ionization modeP^ reconcile fast CM ( < 1 
ps) with the slower band-edge Auger multiexciton relax- 
ation (~ 100 ps) but make the assumption of constant 
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Coulomb matrix elements between X and BX subspaces. 
While this "random- k" -like approximation appears jus- 
tified for the specific case of bulk Si at high energies^ZI, 
it might not apply to direct-gap semiconductor NCs, es- 
pecially since the NC surface is thought to play a large 
role in Auger-like and other Coulomb processes^. Other 
researchers have pursued a generalized treatment, going 
beyond the perturbative approach of typical impact ion- 
ization calculations to examine the effect of phas e and 
population relaxation rates on the CM processes! Cal- 
culations on PbSe NCs using a simple model of the elec- 
tronic structure have shown efficient CM^but, in doing 
so, obtain large values of Coulomb matrix elements that 
appear to conflict with the much slower observed Auger 
rates. The third approach, virtual-state-mediated CM 11 , 
has so far not taken dephasing or population decay into 
account but could be generalized with the introduction 
of decay rates in the energy denominators of their sec- 
ond order perturbation expression (E — ► E — iT). All 
theoretical approaches so far thus hinge on the relative 
rates of coulomb interaction and intraband relaxation at 
threshold, but neither has been measured or accurately 
calculated for NCs. 

In summary, we have determined CM efficiencies in 
CdSe and CdTe NCs by transient photoluminescence. 
Exciton and biexciton features were first characterized 
under 3.1 eV excitation from which we find a relatively 
fast BX radiative rate > 3/c™ d . Measurements un- 
der weak 5.6 and 5.9 eV excitation show no evidence of 
biexciton generation, and thus no CM, up to photon en- 
ergies as high as 3.1E g . 

APPENDIX A: EXPERIMENTAL METHODS 

1. Sample preparation 

Semiconductor nanocrystals (NCs) were synthesized 
by high temperature pyrolysis of precursors^. The prod- 
uct was purified by precipitation once with butanol and 
methanol and redispersed in hexane. Samples were di- 
luted to an appropriate concentration and introduced in 
1mm path length quartz cuvettes. For excitation-flucncc 
studies, dilute solutions with an optical density (OD) of 
w 0.1 at 3.1eV were used, but for measurements under 
UV we employed very concentrated samples (OD w 1 at 
3.1eV) to minimize the effects of degradation. A small 
magnetic bar was added to the solution. It was posi- 
tioned near the excitation point and rotated vigorously 
in a plane parallel to the cuvette surface using an exter- 
nal magnetic stir plate. All measurements and manipu- 
lations were carried out at room temperature. 

2. tPL apparatus and measurement methodology 

The experimental setup used for studying CM by tPL 
with 5.6eV excitation is shown in Fig(5j We isolated the 
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FIG. 5: Apparatus used for studying tPL of NCs with exci- 
tation energies up to 5.6eV 



5.6 eV sum frequency generated beam from the 2.5 eV 
and 3.1 eV fundamental beams using a prism pair. This 
was confirmed by checking that UV-excited tPL signals 
vanish if the relative delay of the fundamental beams 
is adjusted or if a UV absorber (such as a thin glass 
coverslip) is placed in the beam path. For excitation at 
3.1 eV, the OPA was blocked and the 3.1 eV beam was 
directed onto the same excitation spot on the sample. 
For excitation at 5.9 eV, the Ti:sapph was tuned to the 
red and its second harmonic was directly doubled in the 
same nonlinear crystal. The resulting SHG was weak but 
sufficiently strong for tPL collection. 

The excitation beams were characterized spectrally 
with a fiber spectrometer (Ocean Optics QE65000) and 
their spatial profiles were obtained by directly imaging 
onto a CCD camera (Roper Scientific MicroMax) placed 
at the sample position. 

CM determination measurements consisted of alternat- 
ing acquisitions under weak 3.1 eV and UV excitation. 
Exposure to UV was minimized by using low excitation 
power (< 50/iW) and short integration times (~ 5 min). 
Under these conditions, tPL decays excited with weak 3.1 
eV remained unchanged throughout the length of the ex- 
periment, confirming that sample integrity was not com- 
promised. 
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APPENDIX B: DYNAMICS OF MULTIEXCITON 

PL 

1. Kinetic model 

We analyzed the experimentally observed PL dynam- 
ics with a first order kinetic model of multiexciton relax- 
ation: 



dpi (t) 

dt 
dp2(t) 

dt 



dt 



-hPi(t) + k 2 p 2 (t) 
-k 2 p 2 (t) + k 3 p 3 (t) 

-k jPj (t) + k J+lPj+1 (t) 



(Bl) 



where Pj(t) are the relative populations of NCs in the 
j-th multiexcitonic state (i.e. j electrons and holes), and 
kj are the state decay rates. The transient photolumi- 
nescence (tPL) signal seen on the detector will have con- 
tributions from all populated states. We are particularly 
interested in the dynamics of emission from the 15 e — ISh 
transition (i.e. band-edge luminescence), which is de- 
noted by s(t) and given by: 



s(t) ex k\ ad p x {t) + k r 2 ad p 2 (t) + ... = £ k 7 d Pi(t) (B2) 



where k r °" d is the rate of radiation of the j-th multiexciton 
state from the 15-15 level, which in the case of j = 1 and 
j = 2 corresponds to the total radiative decay rates of 
the X and BX, jfc$? d and k r £$, respectively. 



2. Solution of the rate equations 

For clarity, first we solve the model assuming an initial 
population of only X, BX and TX states (i.e. pj(t = 0) = 
for j > 3). We find: 



Ps{t) = P 3 e 
P2(t) = (p° 2 

Pi(t) = (p°x 



-k 3 t 



k 3 - k 2 
k 2 



3 o \ -k 2 t 

Pa e " 



k 



k 2 — ki P2 ' k 2 — ki k 3 - fci 

k 2 ( o , k 3 



h - k 2 
k 2 k 3 



3 -k 3 t 

e 3 



p° e' k ^ 



k 2 - fci V^ 2 ' k 3 -k 2 
k 2 k 3 



Ps e 



-kit 



(h ~ k 2 )(k 3 - k x ) 



P° 3 e- k3t 



where the {Pj} denote initial populations at t — 0. Sub- 
stituting into Eqn. |B2| and grouping terms according to 
their time dependence, we find the luminescence signal is 



given by 



s(t) cx k\ ad (pi + 



k- 2 



k 2 - ki P2 ' k 2 - k x k 3 - ki 

k 3 
k 3 - h 



+ (...)e- fc3t 



Ps e 



We are interested here only in the components of the 
luminescence with dynamics corresponding to the exciton 
and biexciton decay rates. The above result is further 
simplified using w for j > 1 since the multiexciton 
Auger decay rates are more than an order of magnitude 
larger than X decay rates. The decomposition of s(t) as 
a sum of exponentials is thus given by: 



The 



-kit 



-kit 



S(t) 



-kit 



a x = kl ad (j>l + p° 2 + pl) 
asx = (k r 2 ad - k[ ad ) {pl + pl) x 
k 2 P° 3 



1 



k 3 -k 2 p° 2 +pl 



a x e klt + a B xe~ 



k 2 t _|_ / \ g— k 3 t 



component of s(t) is simply proportional 



to the X radiative rate multiplied by the population of 
NCs that start in an X state or higher, since MX states 
eventually decay to the X state. The BX component has 
three factors. It is firstly proportional to the difference 
of X and BX radiative rates because BX luminescence is 
partially offset by a dip in X luminescence, as it is still in 
the process of being populated. The second factor in the 
expression for asx is proportional again to the sum of 
the initial populations in BX or higher states. The third 
term captures the small delaying effect of cascaded decay 
from the TX. In NCs originally in the TX state, there is a 
short t ~ k^ 1 delay before population of the BX, which 
then decays normally at a rate k 2 . Thus, if the e~ k2t 
dependence of the subsequent BX decay is extrapolated 
back to t = one finds a slightly larger value of clbx 
than would be expected in the absence of the delay. 

Fig. [6] serves to illustrate the contributions to s(t) in 
a typical case where p°, p 2 , p 3 ^ 0, as would be the case 
following pulsed excitation and assuming Poissonian pho- 
ton absorption statistics. The solid areas show the de- 
composition of s(t) according to the state involved in the 
emission s(t) = k{ ad Pl {t) + k r 2 ad p 2 {t) + k r 3 ad p 3 {t). The 
dashed lines show the decomposition into a sum of ex- 



ponentials, s(t) 



axe 
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+ (...)e- 



-k 3 t 



as 



would be obtained from a multiexponential fit to a mea- 
sured s(t). 

The expressions for ax and asx derived above for the 
restricted case p° = Vj > 3 suggest a generalization: 



a x = k r 1 ad (p°+p° + ...) = k r 1 ad p° >Q 
asx = (k r 2 ad - k{ ad )p° 3>1 x 



k 2 p° 



(B3) 
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FIG. 6: Decomposition of the tPL decay s(t). Solid ar- 
eas show emission by X,BX and TX states. In contrast, the 
dashed lines show the decomposition of s(t) into a sum of 
exponentials. s(t) = a x e~ kxt + a B xe" kBxt + (. . .)e~ kTxt 



which is proved for the general case of arbitrary {p^} in 
the section lB"5l 



3. Effect of multiexponential X dynamics 

The analysis so far has relied on the assumption of a 
homogeneous sample of NCs with identical decay rates. 
In fact, experimental tPL data often show a distribution 
of X decay rates. Considering that high quality core/shell 
samples do show very close to single exponential dynam- 
ics, we identify two NC subpopulations in each sample, A 
and B, corresponding to high and low quality NCs respec- 
tively. The distribution could be static or mig ht i nvolve 
switching between A and B on slow timescale d 20 | 21 l We 
denote by p A and p^ the relative number of NCs in each 
population in the j-th multiexcitonic state. 

Our experimental data shows that, although the X de- 
cays can be significantly multiexponential, the additional 
contribution attributed to biexciton emission is well de- 
scribed by a single exponential. This suggests that the 
low-quality subpopulation does not contribute to the BX 
signal within our experiment's time resolution. Rapid 
BX decay in those NCs could be explained by the same 
traps or nonradiative mechanisms involved in their fast X 
decay, but could also be due to Auger decay enhancement 
on rough or defective surfaces, as would be anticipated 
from theoretical treatment d < > 29 [ For our purposes we can 
thus approximate the MX decay down to the X state in 
type-£> NCs as instantaneous : 



p*(t > 0) « 



and 



P?(t)=p°* -h(t) (B5) 



where h(t) is an exponential or multiexponential decay 
with h(0) = 1. On the other hand, the A population 



dynamics are single exponential and follow the treatment 
of the previous section. Assuming that the radiative rates 
of A and B NCs are identical, we obtain an expression 
for the tPL signal as follows: 



s(t) cx {k r 2 ad - kl ad )p°* x 1 + 



k-2 



Pj>2 



k 2 p oA 



3>l 



+k[ ad (p^ -e^ t +p^ -h(t))+... 

To simplify, we introduce b — pf ot /p^ t , the ratio of 
the total number of NCs of type B to type A. Since 
photon absorption is not expected to be affected by NC 



quality, P° = bpj. After substituting we discard the 



A superscripts by redefining p^ 
expression: 



S (t) « {k r 2 ad - fer d )pj>i i + 



A _ 



Pj 



and obtain the 



Pi>2 



x k 3 - k 2 p° j>1 
+ kl ad p° J>a (e- k ^ + b-h(t)) + 



(B6) 



where all pj and kj refer to the high quality „4-type NCs. 
Eqn. B6 can be used to interpret the results of our exper- 
imental data analysis. When studying the tPL decays, 
we first summarize the observed X dynamics by fitting to 



a biexponential f(t) = e 



Cfaste' 



Btt , and then 



proceed to fit decays at higher excitation fluence to the 
form s(t) = a,Bxe~ k2t + axf(t) + (■ • ■)■ Comparing to 
Eqn. |B6| we identify 



a BX = (k r 2 ad - fcr d )p?>i i 



Pi>2 



and 



ax = 



1 + 6 

1 + C fast 



h - k 2 P° >x 



y,rad o 
ft l Pj>0 



(B7) 



(B8) 



The first factor in this last expression is sensitive to de- 
tails of sample subpopulations. One expects a priori that 
1, which is the value it would take if X dynam- 



l+b 



l+Cf 



ics were fundamentally biexponential and were captured 



accurately in the fits (i.e. fci = k s i ow and h(t) 



but there could be sample to sample variation to values 
smaller or larger than 1. 



4. Obtaining initial MX populations from tPL fits 

In studying CM in CdSe NCs, we are particularly in- 
terested in determining the fraction of photo-excited NCs 
that started in a BX state. This can be obtained from 
the ratio of fit components agx /a>x i 



7 



obx = k r 2 ad - k[ ad f 1 + b \ p° 3>1 
ax k\ ad \l + c fast J p° >0 



^2 P°j>2 

k 3 - k 2 p° >1 



(B9) 



(a B x /ax) s 



1+A 



r 



where we identify (clex /a<x)sat as the asx/o-x ratio 
one would obtain if only the BX state were initially pop- 
ulated (j>2 = Pj>o)- (cibx I 'ax) sat can be estimated by 
fitting our expressions for asx and ax using values ex- 
tracted from an excitation power series. In an approxi- 
mation, we carried out these fits neglecting the effect of 
the delay term in the asx expression, effectively setting 
A » 0, which should lead to at most a « 25% overesti- 
mate of (asx I 'ax) sat considering the excitation fluences 
used. 

To determine CM yields with tPL, samples are excited 
with weak UV pulses and one looks for a BX signature in 
the resulting decay. Due to energy conservation, p°j >2 = 
for our samples and excitation wavelengths. A simple 
expression is thus obtained for CM yield, y cm , in terms 
of the fit parameters am, «x and the (asx /ax) sat value 
estimated from a previous excitation fluence power series: 



_ P2 



Pl+P2 



asx 
ax 



asx 
ax 



The method described for obtaining y cm is insensitive 
to the precise role of NC subpopulations in X and BX 
decays, and should be valid as long as ax oc p° >0 and 
a BX P°j>i under fairly general conditions. 

On a final note, we see that one can estimate the ratio 
of BX to X radiative rates by: 



h.rad 
K BX 
h,rad 



= 1 



1 + Cf ast \ ( aBX 



aBX 
ax 



Since b is not known independently, the k B a x /k x ld ratio 
thus obtained by assuming b w Cf ast has some uncer- 
tainty. 



C{pj(t)}. Then 

(s + k^P^s) = p° + k j+1 P j+1 (s) (BIO) 
By iterating, one obtains the transformed solution. 
1 



Pn(s) 



s + k n 



Pn 



k n +i o 

Pn+l 



" S + k n+1 

k n +l 



Pn+2 



s + k n+ i S + k n+2 

which has a partial fraction decomposition of the form: 



(Bll) 



Pn(s) 



A 



{n) 



A 



(n) 
n+1 



A 



(n) 



s + k n s + 



Where the A^ are constants independent of s. The 
solution p n {t) are thus given by 

p n (t) = AWe- & »* + A^ 1 e- fc »+ I * + ■ • • = E 4° e_M 



= E — 

o -4- 



j>n 



S + ki 



and the tPL signal is 



j>n 



»(*) = E K^Pnit) = E ( E ^4 



from which we are only interested in the j ' = 1 (X) and 
j = 2 (BX) components: 

s(t) = fc™^! 1 ' e~ M + (fc™ d 4 2) + k\ ad A { ^ e- k2t + ... 



5. Exact result for arbitrary p° 

The system of kinetic equations (Eqn |Bl[ ) can be 
solved for the components with dynamics proportional 
to e~ fc i* and e~ k2t by Laplace transform. Define Pj(s) = 



Expressions for the Aj coefficients involved are found 



by judicious application of the identity 
thijh-jh) in Eqn 



Bll 



(s+a)(s+b) 

In particular, 



k r , 
k n +i - 



J n+1 



k n 

k n +i k r 



r 



k n 
k n +2 - 



Pn+2 



(B12) 



BIO 



that A 



(») 



n+1 



and it can be shown using Eqn. 

1 + r. — fc " ,. I This allows us to reduce the sig- 



nal expression to 
s(t) sa k^A^ e~ klt 



{k r 2 ad 



k[ ad ) 4 2) e~ k2f 
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Where, using Eqn. B12 



_ , O i O i O i 

~ Pi + P 2 + + • • • 



Pj>0 



In simplifying the above two equations we have again 
used the fact that the multiexciton decay rates are much 
faster than the X decay rate. Finally, 



4 2) = P°2 + ( i 



h - k 2 



P.3 



1 4- . \ ^ d • . \ |« 



&3 - &2 



fc 4 - &2 



+ p° + . . .) + (P% + Pi + -..)+• ■ 



P?>1 1 



fc 2 Pj>2 
&3 - fc2P|>i 



which proves the general expressions for ax and olbx 
(Eqn. [B3]fc[B4|. 
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